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ABSTRACT: 



The locations of the fully despun, double synchronous end states of tidal evolution, where 
the rotation rates of both the primary and secondary components in a binary system synchronize 
with the mean motion about the center of mass, are derived for spherical components. For a 
given amount of scaled angular momentum J/ J', the tidal end states are over-plotted on a tidal 
evolution diagram in terms of mass ratio of the system and the component separation (semimajor 
axis in units of primary radii). Fully synchronous orbits may not exist for every combination of 
mass ratio and angular momentum; for example, equal-mass binary systems require J / J' > 0.44. 
When fully synchronous orbits exist for prograde systems, tidal evolution naturally expands the 
orbit to the stable outer synchronous solution. The location of the unstable inner synchronous 
orbit is typically within two primary radii and often within the radius of the primary itself. 
With the exception of nearly equal-mass binaries, binary asteroid systems are in the midst of 
lengthy tidal evolutions, far from their fully synchronous tidal end states. Of those systems 
with unequal-mass components, few have even reached the stability limit that splits the fully 
synchronous orbit curves into unstable inner and stable outer solutions. 

Calculations of material strength based on limiting the tidal evolution time to the age of 
the Solar System indicate that binary asteroids in the main belt with 100-km-scale primary 
components are consistent with being made of monolithic or fractured rock as expected for binaries 
likely formed from sub-catastrophic impacts in the early Solar System. To tidally evolve in their 
dynamical lifetime, near-Earth binaries with km-scale primaries or smaller created via a spin-up 
mechanism must be much weaker mechanically than their main-belt counterparts even if formed 
in the main belt prior to injection into the near-Earth region. Small main-belt binaries, those 
having primary components less than 10 km in diameter, could bridge the gap between the large 
main-belt binaries and the near-Earth binaries, as, depending on the age of the systems, small 
main-belt binaries could either be as strong as the large main-belt binaries or as weak as the 
near-Earth binaries. The inherent uncertainty in the age of a binary system is the leading source 
of error in calculation of material properties, capable of affecting the product of rigidity /x and 
tidal dissipation function Q by orders of magnitude. Several other issues affecting the calculation 
of ^jlQ are considered, though these typically affect the calculation by no more than a factor of 
two. We also find indirect evidence within all three groups of binary asteroids that the semimajor 
axis of the mutual orbit in a binary system may evolve via another mechanism (or mechanisms) 
in addition to tides with the binary YORP effect being a likely candidate. 

Keywords: Asteroids - Satellites of Asteroids - Tides, solid body - Asteroids, dynamics - 
Near-Earth objects 
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1. Introduction 

Over 20 years ago, IWeidenschilling et alJ (|l989l l asked, "Do asteroids have satellites?" Today, 
binary systems have been discovered in every dynamical class of small Solar System bodies 
from near-Earth asteroids to Mars-crossers and main-belt asteroids, among the Jupiter Trojans, 
and in the outer Solar System in the Kuiper belt. Beginning with the Galileo spacecraf t's 



serendipitous disco very in 1993 of tiny satellite Dactyl orbiting (243) Ida (j Chapman et alJ . Il995 



Belton et all . 1 19961 ) while on its cruise to Jupiter and continuing with the success of radar 



lightcurve photometr y, ground-based ada p tive-optics imaging, and Hub ble Sp ace Telescope 
imaging, reviewed by Merline et al. ( 20021 ). Richardson and Walsh (|2006l ). and INoII et aP pOOSi ). 
over 180 small Solar System bodies are suspected to be binary or multiple systems0 In this work, 
we build upon the ideas presented in IWeidenschilling et al. (| 19891 ) to illustrate the tidal end states 
of binary asteroid systems and discuss the inherent material strength of the bodies implied by 
tidal evolution. 

The Keplerian orbit of the two components in a binary system about the center of mass 
allows one to determine the mass of the system and, with an estimate of the component sizes, 
the (assumed equal) densities of the bodies. Density es timates con i bined with taxonomic 



classifications hint at the internal structure of asteroids ( Britt et al. . 20021 ). When compared 



to analog meteorites with known bulk densities, low densities found for rocky asteroids imply 
that, rather than monolithic bodies, asteroids are aggregate bodies made up of many fragments 
with varying degrees of porosity. Aggregates can range from low-porosity fractured or shattered 
bodies to porous, cohesionless rubble piles with "gravitational aggregate" acting as a catch-all 
term for bodies compris ed of many fraginents, i ndependent of porosity, that are held together by 
their collective gravity kicharHson et abl .m. The tidal evolution of these binary systems is 
intimately tied to the internal structure and material strength of the bodies involved, and we will 
exploit this dependence to estimate the combined effect of rigidity and energy dissipation in the 
bodies to determine if asteroids tend to have solid, fractured, or shattered interiors. 

Hints of the internal structure of asteroids also come from the most probable binary formation 
mechanisms for each population: near-Earth binaries likely forr n through rotational disruption y ia 



spin- up by thermal torques due to s unlight (the YORP effect: iBottke et al.l. l2002l: IWalsh et al. 



20081 ) or close planetary encounters (IRichardson et al.l . Il998l : IWalsh and Richardsonl. 120061') . 1- to 



10-km -scale binaries in the main belt may also form through YORP spin-up ( Pravec and Harrii . 
20071). lOO-k r n-sca le main-belt binaries likely form through sub-catastrophic impacts {e.g., 



Durda et al.l. |2004|). and Kuipe r-belt binaries may have formed very early on via gravitational 



collapse (Nesvornv et al. . 20 id ) or later on via collisions or a flavor of dynamical capture 



see 



Nolletal.l .l200d.forarevl;;:i: Binary formation through rotational disr;ption at the very least 



implies a fractured interior, and if satellit e formation i s due to the accretion of smaller particles 
spun off the parent body, as modeled by Walsh et al. ( 20081 ). part or all of the parent body is 



The value of more than 180 suspected binary or multiple systems is taken from 
|http: / / www.johnstonsarchive.net / astro / asteroidmoons.html and based on references therein. Lists 
of binary and multiple asteroid parameters are available from the Ondrejov Asteroid Photom- 
etry Project (http://www.asu.cas.cz/~asteroid/binastdata.htm) and the Planetary Data System 
(|http://sbn.^sredu/pds/asteroid/EAR_A_COMPIL_5-BINMP_V3-0/data/D. 
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likely an aggregate of smaller pieces and structurally weak. For sub-catastrophic impacts, the 
satellite is either a solid shard or an aggregate made up of impact ejecta while the parent body 
may remain solid or fractured rather than shattered. Thus, we anticipate that examination of 
tidal evolution in binary asteroid systems will find that large main-belt binaries are consistent 
with solid or fractured rock, while near-Earth and small main-belt binaries are consistent with 
structurally weaker, shattered or porous gravitational aggregates. 

Section [2] introduces the important dynamical quantities of binary systems, angular 
momentum and energy. In Section [3l we derive and illustrate for a range of angular momenta the 
fully despun, double synchronous orbits that are end states of tidal evolution and then discuss 
stability limits and energy regimes in Section HI Using the known properties of binary asteroids in 
the near-Earth region and the main belt, in Section [5] we constrain the material properties of the 
bodies based on tidal evolution. 



2. Angular momentum and energy content 

The dynamics of a binary system with primary of mass Mp and secondary of mass Mg in 
a mutual orbit about their common center of mass can be described equivalently by a system 
where a body of mass MpMs/(Mp + Mg) orbits a stationary mass, Mp + Mg. The orbital angular 
momentum of such a system is L = AfpMs/(Mp + Mg) \J G (Mp + Mg) a (1 — e^) where G is the 
gravitational constant and a and e are the semimajor axis and eccentricity of the mutual orbit, 
respectively. Defining the component mass ratio q = Mg/M-p = {p^/ Pp){Rs/Rp)^ < 1, where p 
and R are the bulk density and radius, respectively, of the components, using the system mass 
Mgys = Mp + Ms = (l+q) Mp, and applying Kepler's Third Law, n^a^ = G (Mp + Mg) in terms 
of the mean motion n, the orbital angular momentum is: 

L = ^^Mg,ga'n{l-eY' = Mp a' n {l - e^' ■ (1) 

For principal axis rotation, the spin angular momentum of the system is given by the moments of 
inertia / = aMR'^ of the components and their spin rates w as: 

+ Igo^g = a,M,Rlu;Jl + ^(P^Y\'/'^ 



S = IpUp + hcos = apMpR'pOJp\^l + ^i^^j q'^' ^ j ■ (2) 

The coefficient a is 2/5 for a uniform-density, rotating sphere, but varies with the shape of the 
body and the density profile of the interior. 

For comparison, suppose a rapidly spinning, uniform-density parent body sheds mass in such 
a way to conserve angular momentum and produce the aforementioned mutually orbiting binary 
system. When a spherical body with mass M and radius R spins at the breakup rate Wbreak 
without cohesion among its constituent particles, the inward acceleration due to gravity GM/R? 
at the equator is balanced by the outward centrifugal acceleration i^brcak-^ rotation such 

that Wbreak = yGM/R? . The angular momentum J contained in the critically rotating sphere is 

/Wbreak Or: 



J = - ^GMg3ygi?eff, (3) 
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where the mass and radius of the parent body have been written as the total mass and effective 



1/3 



radius of the binary system, R^g = {R^ + Rg)^^^ = yl + ^ qj Rp, produced by the breakup 
of the parent body. Subsequently, the total angular momentum of a binary system J = L + S 

IgMI R^Q (attributed to IParwinl (|l887l )) such that J / J' ~ 0.4 



is often normalized by J' - . , ^^.^^yg^ 

indicates the binary could have formed by ni a ,ss sh edding from the spin-up of a single spherical 
strengthless parent body. iPravec and HarrisI (j2007l ) normalize the total angular momentum of 
the binary system by that of the critically rotating spherical parent body (Eq. [3]) such that their 
scaling parameter = (J/J')/0.4, and ol = 1 has the same implication for binary formation. 
We utilize the J' normalization throughout this work. 

The minimum and maximum separations of two components in a binary system are limited by 
the physical size of the components and the total angular momentum of the system, respectively. 
The separation of two components is naturally bounded from below by the contact condition 
where, at a separation of i?p + i?s, the components are resting against one another. In terms of 
the semimajor axis and the radius of the primary, the contact limit is: 



a 

lir 



1 + 



R^ 

Rr. 



1 + 



Pp 



1/3 



(4) 



^P/min ^"-P 

The contact limit ranges from a/Rp = 1 to 2, akin to a pea resting on the surface of a basketball 
to two basketballs in contact. One could consider the Roche limit as a soft lower bound on the 
separation of the bodies, as opposed to the hard limit of physical contact, but with a modest 
amount of cohesion, a secondary may exis t with i n the Roche limit and dow n to the contact limit 
without disrupting (iHolsapple and Michel bnnsi : iTavlor and Margotl boid ) . An upper bound is 
placed on the separation of the components by the total angular momentum content of the system. 

If the entire budget of angular momentum J = (yj / J'^ ^ GM^y^Rcs is transferred to the circular 
mutual orbit of the components, which are no longer spinning (S" — t- in Eq. [2]), the maximum 
attainable separation according to Eq. [T]with e = is: 



a 



{J/J' 



n2 {l + qy 



1/3 



(5) 



•-P / max y \ Ps 

which increases quickly with decreasing q, allowing for binaries with smaller secondaries to have 
much wider separations than binaries with similar-size components for the same amount of angular 
momentum. 

Ignoring external influences on the system, the total energy of the binary E is comprised 
of the rotations of the components (1/2) /w^, the orbital motion about the center of mass 
(1/2) MpMs/ (Mp -|- Ms) a^n^, and the mutual gravitational potential for a two-sphere system 
-GMpMs/a such that: 



E 



Wp + 



2 



1 MpMs 
2Mp + Ms 



2 2 

■an. 



(6) 



Whil e shape deformation may occur due to changes in spin rate ( Harris et al. . 20091 : Holsappl3 . 
20inl ) as the binary system evolves, we assume that this does not result in significant interchange 
of self-potential energy with the terms in Eq. [6l The energy accounted for in Eq. [6] is free to be 
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traded among its components via the dynamical interaction of the two bodies, and the total energy 
can be dissipated as heat as a result of internal friction due to tidal flexure. Depending on the 
amount of rotational and translational energy of the system compared to the mutual gravitational 
potential energy, the total energy may either be positive or negative. A negative total energy 
requires the system to remain bound, while a positive total energy could allow the components to 
unbind themselves through their dynamical interaction. 



3. Fully synchronous orbits 

The natural evolution of a binary system is through mutual tidal interaction. The 
differential gravity across each body due to the proximity of its companion acts to evolve 
the rotation states of the individual bodies as well as their physical separation; tides tend 
to evolve rapidly rotating components in close proximity to a more widely separated system 
with more slowly rotating components. The end state of tidal evolution for a binary system 
is a fully despun, double synchronous orbit (hereafter called a fully synchronous orbit), where 
the spin rates of both components, Up and ujg, equal the mean motion in the mutual orbit n. 
The classic example of such an end state is in the Pluto-Charon system where the bodies are 
face- locked ( Christv and Harrington . 1978l l. meaning they keep the same sides facing one another 



because the periods of rotation and revolution have synchronized. 

The scaled angular momentum of the system J /J' at any time is given by the sum of Eq. [1] 
and Eq. [2] divided by J'. Assuming a circular mutual orbit, which the majority of binary asteroid 
systems have, and after setting the spin rates ujp and ujg equal to the mean motion n, replacing 
n with a function of the semimajor axis a via Kepler's Third LawH and some rearranging, the 
locations of the fully synchronous orbits agync/Rp are the solutions to the quasi-quadratic equation: 



Oip \Pb J ^ Ctp V Ps / 

If the components are spherical and have similar uniform densities, the condition for a fully 
synchronous orbit reduces to: 

5 q 1 fa^y^^Y + + 1 = 0. (8) 



2 1 + g 1 + g5/3 \ Rp J 2 1 + g5/3 v / ^ \ ^ 

Depending on the total angular momentum of the system J/ J' and the mass ratio q, the equation 
above may have two, one, or zero real solutions corresponding to inner and outer fully synchronous 
orbits, a single degenerate fully synchronous orbit, or the absence of a valid fully synchronous orbit, 
respectively, for the system. Solutions to Eq. [8] for all q are shown in Fig. [1] for various values of 
J/ J' where the number of solutions for a specific q is given by the number of intersections between 
the curve of the relevant J/ J' value for the system and the horizontal line of constant mass ratio 



^In using n?a^ = G {Mp + Mg) here, we imply the components of the binary system have 
spherical shapes or may be treated as point masses. 
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q that the system will tidally evolve along. From Fig. [H binary systems with smaller secondaries 
can clearly attain much wider separations than binary systems with components of similar size 
as predicted from the angular momentum limit in Eq. [5j Larger amounts of angular momentum 
are required to support systems with larger secondaries, and as angular momentum increases, the 
synchronous orbit curves surround one another with the inner synchronous orbit moving inward 
while the outer synchronous orbit moves outward approaching the angular momentum limit for 
smaller mass ratios. 

Figure [J is a hybrid of the "universa l diagram" for tidal evolution of ICounselmanI (jl973l ) 
and Fig. 1 of Weidenschilling et al. ( 19891 ) in that it combines the ability to directly read off the 
fully synchronous orbits for specific values of angular momentum of Counselman with the mass 
ratio dependence of Weidenschilling et al. Furthermore, our figure accounts for the spin of the 
secondary, which Counselman does not, and can be modified to include tidal evolution timescales 
as done by Weidenschilling et al., which we will do in Section [5l Understanding evolution through 
our version of the universal diagram is aided by Fig. [2l which illustrates the general trends of 
tidal evolution when the spins of the components have the same sense of rotation as the motion 
in the mutual orbit, i.e., all motion is prograde. If the configuration of the system falls between 
the two solutions to the fully synchronous orbit equation for the q of the system, meaning it 
visually appears "below" the fully synchronous orbit curve in ((7, a/i?p)-space like the lower arrow 
in Fig. [21 the system will evolve outward as spin angular momentum is transferred to the mutual 
orbit via the tidal interaction (wp > n, ~ n). If the system sits "above" its fully synchronous 
orbit curve in the position of one of the upper arrows in Fig. [21 due to having a large secondary 
or a wide separation, the orbital angular momentum makes up a larger fraction of the total 
angular momentum such that the angular momentum available in the spins of the components 
requires that Up < n (for ~ n). The system must then evolve inward as angular momentum 
from the orbit is transferred to the spins of the components. For systems with constant angular 
momentum, the cases of inward evolution require a binary formation mechanism that initially 
produces well-separated components with rotation rates slower than the mean motion because 
tidal evolution cannot have evolved a system outward to these configurations above the fully 
synchronous orbit curve. Note that binaries with equal-mass {q = 1) components in Fig. [Tjonly 
have fully synchronous end states if J/J' > 0.44. For J/ J' < 0.44, cases can exist where q lies 
entirely above the fully synchronous orbit curve (the uppermost arrow in Fig. [2]) meaning the 
fully synchronous orbit equation has no solution for that value of q. The resul ting inward ti dal 
evolution may be important in the context of forming contact binary asteroids ( Taylor . 20091 ) and 
will be addressed in the future. 



4. Stability and E = limits of fully synchronous orbits 

The fully synchronous orbit solutions to Eq. [7] are equivalently thought of as contours of 
constant J/ J' in {q, a/i?p)-space satisfying: 



\ 2/3 



T 1 / \ 1/2 l + q5/3 

J _ 1 q I ^ \ o!p \ps J ^ 



+ «p . '^V. ^ • (9) 



a 



-3/2 



The maxima of the set of contour curves for a continuous range of J/ J' are traced out by taking 



7 



2 4 6 8 

Separation, a/Rp 



Fig. 1. Component separations async/Rp for the fully synchronous orbits of binary systems of 
given a mass ratio q and scaled angular momentum J / J' . Solid lines indicate the solutions to 
Eq. El for J /J' = 0.25, 0.325, 0.4, 0.425, 0.44, 0.45, 0.475, and 0.5. The shaded regions indicate 
separations that are inaccessible to the binary components. The region at left (Eq. U]) requires 
the components to be in contact, while the region at right (Eq. [Sj shown only for J/ J' = 0.5) is 
disallowed by angular momentum conservation for circular orbits. The dashed line indicates the 
synchronous stability limit (see Section U]) that divides the two solutions for each mass ratio, if 
they exist, into an unstable inner synchronous orbit and a stable outer synchronous orbit. In most 
cases, the unstable inner synchronous orbit is within the contact limit. The dotted line indicates 
where fully synchronous systems have total energy E = (also see Section [4]). 
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Fig. 2. Directions of tidal evolution for systems with J / J' = 0.4. For binaries dominated by the 
spin angular momentum of the primary mass, systems under the fully synchronous orbit curve 
will evolve outward as angular momentum is transferred from the rapid spin of the primary to 
the orbit (cjp > n). Systems above the fully synchronous orbit curve evolve inward as angular 
momentum is transferred from the orbit to the spins of the components {uj^ < n) until either the 
outer synchronous orbit is reached or the orbit collapses to contact. The shaded regions indicate 
separations that are inaccessible to the binary components due to the contact limit and angular 
momentum limit for J/J' = 0.4. The dashed line indicates the synchronous stability limit, while 
the dotted line indicates where E = for fully synchronous orbits (see Section [4]). 
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the derivative of Eq. [9] with respect to a/Rp giving: 



Ostab 
Rn 



3ar 



1 + 



1 + 



Or 



2/3 



5/3 



PsJ 



1/2 



For spherical components with similar uniform densities: 



Q^stab 
Rn 



6 i_ 

5 



5/3 



.1/2 



(10) 



(11) 



which is the dashed curve in Fig. [T] that splits the sy nchronous orbit curves into inner and outer 
solutions when there are two real solutions to Eq. [8l iHarris and WardI (119821 ) showed that Eq. [10] 
is also a stability limit against perturbations to the system such that the outer synchronous orbit 
in the two-sphere case that satisfies: 



a 

Rr, 



> 



6 1 + 
5 ^ 



1 



5/3 



1/2 



(12) 



is stable, while the inner synchronous orbit is not, meaning that the binary system will tend 
to evolve away from the inner synchronous orbit given the opportunity. If perturbed outward 
from the inner synchronous orbit, the components will tidally evolve outward to the stable outer 
synchronous orbit (the lower arrow in Fig. [2|); if perturbed inward, the components will fall to 
contact, though for most J/J', the inner synchronous orbit is already within the contact limit. 
The arrows in Fig. [2] also show that perturbation of the system from the outer synchronous orbit 
will allow the system to return to the outer synchronous orbit, indicating its stability against 
tidal perturbation. It is also clear from Figs. [1] and [2l that the stability limit is always beyond 
the contact limit such that whenever a particle is lofted from the surface of the primary or when 
a parent body fiss ions, the compo nents are initially in a tidally unstable configuration that may 
lead to re-impact (IScheeresl - liooi l. 



For bodies with «p,s 7^ 2/5 , one may use the parameter Ap^s = 5ap^s/2 introduced 



by Descamps and MarchisI ( 20081 ) to find that the stable outer synchronous orbit is described by: 



a 

Rn 



> 



6 1 + g 
5 a 



An + As 



,5/3 



1/2 



(13) 



correcting both the formula and direction of the inequality given by iMarchis et al. (l2008al lbl). We 
note, though, that if one (or both) of the components of the binary system is nonspherical, Eq. [13] 
is not completely accurate because the formulation of the fully synchronous orbit curves in Eq. [7] 
and thus the stabilit y limit, would b e complicated by the mutual gravitational potential between 
nonspherical bo dies (IScheerea. 20091) a nd the synchronous rotation rate required for relative 
equilibrium (see Scheeres, 2006 . 2007al ) . Therefore, Eq. [13] is exact only in the two-sphere case 



with Op^s 7^ 2/5 due to density inhomogeneity within the individual components. 

The stability limit in Eq. [10] gives the minimum separation of a stable, fully synchronous 
tidal end state for a given mass ratio q. The corresponding minimum angular momentum required 
for a system to have a stable tidal end state results from substituting Eq. [10] into Eq. [9] and, by 
assuming similar component densities and Op^s = 2/5, can be written as: 
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VI35J (l + <;)23/12 



l + q 



5/3 



1/4 



(14) 



This critical value for the angular momentum is equivalent to that derived by iHut m 
the context of binary star systems. Satisfying the stability condition is a necessary, but not 
sufficient, condition for a stable end to tidal evolution. Systems that do not satisfy the stability 
condition have not reached a stable tidal end state, while systems that do satisfy the stability 
condition may or may not continue their evolution depending on the angular momentum of the 
system. At the minimum angular momentum of Eq. [T^ the synchronous orbit equation (Eq. [8|) 
has one degenerate solution, which is simply the stability limit, and tidal evolution may cease. 
If the angular momentum of the system exceeds J/ J^^.^^. , the synchronous orbit equation has two 
solutions, and the system will continue to evolve to the outer synchronous orbit. Of course, 
having less than J / J^j,j^ will cause the system to collapse as the synchronous orbit equation has no 
sol ution regardless of w hether or not the stability condition is satisfied. This argument is used 
bv lLevrard et al.l (|2009l ) to posit that nearly all transiting extrasolar planets, in the absence of a 
parking mechanism, will inevitably collide with their host stars because these planetary systems 
lack sufficient angular momentum to support stable tidal end states for their components. 



Upon full synchronization (wp 



n] 



the total energy of the system (Eq. [6]) may be 



positive or negative with the E = limit occurring at: 



OE 
Rr, 



1 + 



1 + 



Os I' f>p\ 
Ps) 



a 



2/3 



5/3 



1/2 



(15) 



where spherical shapes have been assumed in determining the mutual gravitational potential. 
Fully synchronous orbits with semimajor axes Ogync > OE have E < because the combination 
of translational and potential energy (oc falls off more slowly with increasing a than the 
rotational energy of the components (oc a^^^ when synchronized. Fully synchronous orbits with 
ctsync < CLE h ave -E > and could evolve to escape given sufficient perturbation from relative 
equilibrium ( Bellerose and Scheeresl . boosi ). Note that in the two-sphere case oe = Ostab/v^ 



such that all stable, fully synchronous orbits have £^ < 0, but the converse, having E < 0, does 
not necessarily guarantee stability. For systems with J/J' > 0.334 in Fig. [H the realistic fully 
synchronous orbit solutions, i.e., those beyond the contact limit, have E < 0. Systems with less 
angular momentum can have unstable inner synchronous orbits for which E > 0. 

As an interesting side note, if one inserts ce into L and S in Eq. [1] and Eq. [2] for a fully 
synchronous orbit, one finds that L/S = 1 or the angular momentum of the system re volving 
abou t its center of mass is equal to the rotational angular momentum of the components (jScheeresl . 
20091 ). Similarly, inserting Ostab finds L/S = 3 01 t he angular momentum of t he system's revolution 
is three times the rotational angular momentum ( Hut , 198d : Scheeres . 20091 ) . Thus, a stable fully 
synchronous tidal end state must satisfy L/S > 3 or, equivalently, the orbital angular momentum 
must account for more than 3 /4 of the total angula r momentum of the systeni, which is obser ved 



in ne arly equal-mass binar y systems (90) Antiope (jMichalowski et al.l . l2004l : iDescamps et al. 



2OO7I ) and (69230) Hermes dMargot et al.[ bood ) 
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5. Material properties of binary asteroids 

During the tidal evolution process, angular momentum is transferred between the spins of the 
components cOp and and the mutual orbit abo ut the center of mass of the s ystem. For a progra de 
syste m, angular momentum conservation {e.g., Murrav and Dermott . 19991 : Taylor and Margot . 
20ld ) requires that the change with time of the orbital separation of the components follows: 



a 



X 



19 fipQp 
sign (wp 



9(1 + 



a/2 



a 

R„ 



N , Rs MpQp ■ / 
^) + -^ ~n ^^^^ '-^'^ 



-11/2 



n] 



(16) 



The sign of the quantity uj — n for each component determines whether the separation increases 
or decreases due to tides raised on that component; if > n (w < n), tides slow (accelerate) 
the rotation of the component and increase (decrease) the orbital separation. The contribution 
from tides raised on the secondary are a factor of the size ratio weaker than from tides 
raised on the primary and, furthermore, the smaller secondary will tend to synchronize to the 
mean motion n before the primary causing the tidal torque on the secondary, as well as its 
contribution to the change in orbital separation, to vanish earlier. The overall strength of the 
tidal mechanism is inversely related to the product fiQ, which combines the rigidity or shear 
modulus fi of the materials and the tidal dissipation function Q that relates t o the lag angle 
between the axis of sym r netry of the tidal bulge and the tide-raising body (see iGoldreichl . Il963l : 
Efroimsky and Williams! . I2nn9l . for detailed discussions). The larger the quantity //Q, the more 
resistant the material is to orbit evolution due to the tidal mechanism. Here, for completeness, we 
have allowed for the primary and secondary to have their own fiQ values. 



5.1. Caveats 

In estimating the material properties of binary asteroids, we assume that tides are the 
dominant method of dynamical evolution in these systems, specifically in the evolution of the 
separation between the components. It has been argued that the binary YQRP efi'ect (BYORP; 
Cuk and Burns! . 2005 : Cuk and Nesvornv . 2010 : McMahon and Scheeres . 2010a ). where a 
synchronous secondary acts to asymmetrically re-radiate sunlight with respect to its orbital 
velocity so that the orbit i s exp a nded or contracted , can act on timescales faster than tid al 
evolution (|Cuk and Burns! . 1200,4 lOoldreich and S ari'. '2009"; 'Mc Mahon and Scheere^ . l2010bl ). 
This mechanism is similar to the YORP effect (jRubincam, 200Q), an asymmetric re-radiation 



of su r ilight with respect to rotational velocity. Ho wever, unlike the YORP effect (iTaylor et al 



20071 : iLowrv et all . l2007l : iKaasalainen et al.l . ^OOT! ) , BYORP has yet to be proven observationally. 
Furthermore, BYORP requires a synchronous secondary and thus cannot be the dominant active 
evolution mechanism for systems with asynchronous secondaries. Also of interest, specifically 



^In deriving Eq. [T6| it is assumed that for both components the rigidity fi dominates the stress 
due to self-gravity pgR ^ Gp^R^, where g is the surface gravity, which is reasonable for bodies 
under 200 km in radius ( Weidenschilling et al. . 1989! ) such that the potential Love number k2 oc p~^. 
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among small binary astero ids with primaries less than 10 km in diameter, is the idea of mass 
lofting or "tidal saltation" ( Harris et al. . 20091 ). where particles at the equator of a rapidly spinning 
primary become weightless due to the gravitational presence of the secondary passing overhead, 
briefly enter orbit, and transfer angular momentum to the orbit of the secondary before falling 
back to the surface of the primary. This method has also be e n arg ued to expand the mutual orbit 
more quickly than tidal evolution (jFahnestock and Scheeresi. l2009|). It may be a combin ation of 
effects, also including close planetary flybys (Farinella, 1992; Bo ttke and Melosh . 1996al jbl). that 
evolve the separation of near-Earth binary components, but tides are the only mechanism one 
can say at this time must act on all systems. For main-belt binaries with 100-km-scale primaries, 
tides should be the dominant mechanism because close planetary encounters are not feasible in 
the main belt, the thermal YORP and BYORP effects are weakened by the increased heliocentric 
distance to the main belt (and the increased sizes of the bodies involved), and the primaries do 
not rotate rapidly enough nor are the secondaries close enough to produce mass lofting. 



5.2. Estimation of fiQ 

Integration of Eq. [16] provides, symbolically, the evolution of the orbital separation as a 
function of time in terms of the product /ipQp and the time over which tidal evolution has taken 
place At. Both f^pQp and At are inherent unknowns in any binary system, but their ratio is fully 
determined by the current separation of the binary system a/Rp, the physical properties of the 
components, and the assumption of an initial separation of the components after the formation 
of the binary. The ratio is very weakly dependent on the initial separation (see Section [5.6p : we 
choose 2Rp because it is both the contact l imit for equa l -size components and a reasonable initial 
separation for binaries formed via spin- up ( Walsh et al. . 20081 ) . Thus, the evolution of the orbital 
separation gives: 



fJ-pQp _ 52^/3 



At 



19 



a/2 



a 

Rn 



13/2 



2I3/2 



^ ^ Rs i^'pQp ''"s I 
Rp fisQs At J 



[17) 



for a system where tides on both components are causing the mutual orbit to expand with Tg 
representing the length of time the secondary contributes to the tidal evolution, which is less than 
or equal to the age of the binary system At. Ignoring the rightmost term in parentheses for the 
moment, to estimate either fipQp or At, the other must be assumed, where varying one by, say, an 
order of magnitude varies the other parameter by the same amount. With a judicious assumption 
of the age of the binary , one can estimate UpQp for the primary or at least place rough bounds on 
its value ( Margot et al. . 20021 : Mar got and Brown . 2003); simila r ly, by assuming a va l ue of U pQp, 
one can estimate the age of the binary (jWalsh and Richardsonl . l2006l : iMarchis et al.l . l2008allbl: 
Goldreich and Sari . I2OO9I ). A related analysis was done for the martian moon Phobos by lYoderl 
(|l982l ). who estimated its fiQ is at least of order 10^^ Nm ^ based on the eccentricity evolution of 
Phobos within the At inferred from the minimum age of its surface. 

If one ignores the contribution of the secondary, equivalent to setting Tg to zero, Eq. [T7] 
represents a lower bound on fipQp for a given timescale At. Allowing the secondary to contribute 



13 



for the entire age of the binary by letting Tg — ?■ At gives an upper bounc^ on /ipQp, which is a 
factor of 1 + (Rg/Rp) {fJ-pQp/ fJ-sQs) greater than the lower bound found using only the primary. 
If At is correct, the actual fipQp should lie between the two bounds as the secondary tends 
to be synchronously rotating such that Tg < At rather than contributing to the mutual orbit 
expansion over the entire age of the system. Ideally, one would know when the secondary became 
synchronized, evolve the system using both tides until Tg, then include only the tide on the 
primary for the remainder of At. Our ignorance of if/when the secondary became synchronous 
cannot change the estimate of HpQp by more than a factor of 1 + (Rs/Rp) {iJ-pQp/ fJ-sQs), which is 
nominally less than two for components with equivalent material properties, far less critical than 
the choice of At, the age of the binary, but similar to the effect of measurement errors on the 
physical parameters in Eq. [17] (see Section 15. 6p . 

In examining orbit expansion by tides, one cannot directly determine whether the primary and 
secondary have different nQ values, though having a stronger or weaker secondary can change the 
estimate of fJ-pQp- Assuming equivalent nQ values for the components, fipQp cannot be affected 
by more than a factor of two by accounting for the secondary; often the factor will be much less 
due to the size ratio of the components and the fact that synchronous secondaries require Tg < At. 
If the secondary is mechanically stronger than the primary such that UsQs > IJ-pQp, the effect on 
HpQp compared to ignoring the secondary in Eq. [T7] must still be less than a factor of two because 
(Rg/Rp) ifJ-pQp/ fJ-sQs) (Ts/^t) < 1. Thus, having a stronger secondary than the primary serves to 
reduce the upper bound on fipQp for a specific age At, but the actual material strength of the 
secondary is not constrained. On the other hand, if the secondary is mechanically weaker than the 
primary such that UsQs < ^^pQp^ depending on the size and ratios of the components and the 
age of the system At, the effect on /ipQp compared to ignoring the secondary could conceivably 
exceed a factor of two, but is unlikely to differ by more than an order of magnitude unless the 
system is very young and the secondary is orders of magnitude weaker than the primaryj^ 

For a synchronous s econdary, the eccentricity can increase or de crease depending on how 
/isQs compares to HpQp (Goldreich, 19631 : Goldreich and Soter . 19661 ) with the prediction that, 
for equivalent material properties , systems with monolithic components and q < 0.31 will have 
the eccentricity excited via tides ( Harris and Ward . 19821 : Weidenschilling et al. . 19891 ). Given we 
observe the vast majority of binaries, including those with much smaller mass ratios than 0.31, 
have circular mutual orbits, to prevent eccentric i ty exc itation the secondaries must be mechanically 
weaker than the primaries ( Mar got and Brown . 20031 ) by roughly an order of magnitude to lower 
the limiting mass ratio from 0.31 to below the observed mass ratios in the binary population, 
e.g., by a factor of 30 for the smallest mass-ratio system we will discuss, (702) Alauda. However, 



^We do not use the terms "lower bound" and "upper bound" in a rigorous mathematical sense 
here because factors such as uncertainties in the physical properties used in Eq. [T7] can affect the 
values of these bounds by factors of order unity. 

^The time it takes for the secondary to synchronize Tg shortens compared to At as ^sQs is made 
weaker so that Tg/At tends to cancel out the effect of fJ-pQp/ fJ-sQs growing. Only if Tg and At are 
comparable due to the youth of the system will the disparity between the mechanical strength of 
the primary and secondary have a strong effect on the calculation at hand (if it can overcome the 
size ratio factor as well). 
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Goldreich and Sari suggest that binary systems with gravitational-aggregate components, 

rather than monohthic ones, do not require much weaker secondaries to have circular mutual 
orbits as the eccentricity will damp for all mass ratios if the components are aggregates of similar 
material. So again, it is difficult to constrain the separate values of /xQ for the components, but our 
ignorance should not strongly affect our results for particular systems unless the extreme case of a 
very young binary with components having very disparate material properties exists. Therefore, 
when estimating the material properties of a binary system via Eq. [T71 we may ignore the 
possibility of the components having different material properties to obtain an order-of-magnitude 
result Ifl In the following sections we will also ignore the contribution of the secondary in Eq. [17] 
to find a lower bound for //Q of the system with a specific age At. An upper bound on /iQ from 
accounting for the secondary is nominally less than a factor of two larger if the material properties 
of the components are similar. 



5.3. Large main-belt binaries 

First we consider binary systems in the main asteroid belt and among the Jupiter Trojans 
with 100-km-scale primary components that have been characterized by direct adaptive-optics 
imaging and lightcurve photometry. To be included in Table [H orbital properties must be known 
and some estimate of the sizes must be available. Most of the binary systems discovered via 
imaging of large main-belt asteroids and Jupiter Trojans have secondaries roughly one-tenth 
the size of their primaries {q ~ 0.001) or smaller. These secondaries are likely the result of a 
su b-catastrophic impa ct on the parent body {i.e., SMATS, smashed target satellites, as described 
by Durda et al. ( 20041 )). which is supported by the angular momentum budget of the systems. 
With the exception of (90) Antiope and (617) Patroclus, which have nearly equal-size components, 
the 100-km scale main-belt binaries have J/ J' values of roughly 0.2, well below the J/ J' ~ 0.4 
regime characteristic of binaries formed via a spin-up mechanism. Among the q < 0.1 binaries, 
the average primary component spins with a period of roughly 6 h, twice the rotational breakup 
period, yet accounts for 97% of the angular momentum of the system due to the large rnass 



disparity between t he components. Nearly equa l-mass binary (90) Antiope (jMerline et al.l . l200d : 
Michalowski et al.l . bood : IPescamps et al.1 . l2007l ) has J/ J' ~ 0.5 that is more similar to a binary 
formed through spin-up than through a sub-catastrophic collision. However, the sheer size of 
the components, each more than 80 km in diameter, and its location in the main belt make it 
difficult for the YORP effect or close planetary encounters to explain how the (90) Antiope system 



2001; Marchis et al 



originally formed. Nearly equal-mass binary (617) Patroclus (jMerline et al 

20061 : iMueller et al.l . l2010i ) has J/ J' ~ 0.8, far larger than the other large main-b elt binaries we 



consider, but near the upper limit for giant impacts of J/ J' < 0.8 ( Canup . 20051 ). Such a high 
angular momentum content is r nore similar to m any of the Kuiper-belt binaries that may have 
formed via n-body capture (see iNoll et al.l . l2008l . for a review) rather than spin-up or collisions. 

The q < 0.1 binaries in the main belt are in the midst of a lengthy tidal evolution. As 
shown in Fig. [3l of the q < 0.1 binaries, only (22) Kalliope has reached the synchronous stability 
limit, but even (22) Kalliope is very far from reaching the outer synchronous orbit that nearly 



^By ignoring any possible difference in fiQ of the components, Eq. [T71 gives the value of fiQ for 
the system, as in it applies to both components rather than the primary only. 
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Name 


Rp (km) 


Rs/Rp 


q 


p (g cm ^) 


a/Rp 


J/J' 


fiQ (Nm-2) 


22 


Kalliope 


85 


0.213 


0.009664 


2.5 


12.5 


0.22 


3.3 X 10^2 


45 


Eugenia 


98 


0.036 


0.000047 


1.1 


12.1 


0.19 


3.4 X 10^ 


87 


Sylvia 


128 


0.063 


0.000250 


1.5 


10.6 


0.20 


1.6 X 10^1 


90 


Antiope" 


43 


0.955 


0.871 


1.3 


3.9 


0.49 


<4.0 X 10^6 


107 


Camilla 


103 


0.050 


0.000125 


1.9 


12.0 


0.18 


4.1 X 10^° 


121 


Hermione 


103 


0.066 


0.000287 


1.1 


7.5 


0.22 


5.1 X 10^^ 


130 


Elektra 


90 


0.026 


0.000018 


3.0 


14.0 


0.14 


5.1 X 10^ 


283 


Emma 


73 


0.079 


0.000493 


0.8 


8.2 


0.20 


1.1 X 10^^ 


617 


Patroclus'' 


51 


0.920 


0.779 


1.3 


13.5 


0.83 


1.5 X 10^3 


702 


Alauda 


97 


0.018 


0.000006 


1.6 


12.6 


0.13 


7.9 X 10® 


762 


Pulcova 


67 


0.160 


0.004096 


1.9 


12.2 


0.17 


5.1 X 10^^ 



Table 1: Physical properties of main-belt and Jupiter- Trojan binary asteroids with primary radii of 
order 100 km. Rp is the radius of the primary; Rg/Rp is the size ratio of the components; q is the 
mass ratio, taken to be the size ratio cubed by assuming components of equal density p; a/Rp is the 
semimajor axis of the mutual orbit; J/ J' is the normalized angular momentum assuming spherical 
bodies and a synchronous secondary using Eqs. [T] and [21 for which an average value of roughly 0.2 
for q < 0.1 binaries indicates a sub-catastrophic collisional origin rather than formation through 
spin-up. For the triple systems (45) Eugenia and (87) Sylvia, only properties of the larger outer 
satellite are considered. Values of pQ correspond to tides raised on the primary acting over the 
4.5 Gyr age of the Solar System. Accounting for tides raised on the secondary would increase pQ 
by no more than a factor of two. Data are taken from the Ondrejov Asteroid Photometry Project 
binary asteroid parameters table (http://www.asu.cas.cz/~asteroid/binastdata. htn?, 2010 April 8 
releas e) with the exception of (702) Alauda whose parameters are provided by iRojo and Mar got 
(|201lh . 

(90) Antiope is known to be in a fully synchronous state where tidal evolution has ceased. Because 
of its nearly equal-size components in such close proximity, its tidal evolution occurred over much 
less than the age of the Solar System. 

^ (617) Patroclus is believed to be in a fully synchronous state and must have formed at a larger 
initial separation than the value of 2Rp assumed here. 
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coincides with the angular momentum hmit for these systems. The inner synchronous orbit for 
J/ J' = 0.2 hes above the contact hmit, so these binaries must have begun their tidal evolution 
from an initial separation of rough ly 2R^ or more to have evolved outward due to tides. Nearly 



equal-mass binaries (90) Antiope (IMichalowski et al.l . 1200J; iDescamps et al.l . 120071 ) and (617) 



Patroclus (jMarchis et al.1 . boOft iMueller et al.1 . l20ld ) are both believed to be in their fully 



synchronous tidal end states where the rotational periods of the components equal the period of 
the mutual orbit. Of the large main-belt binaries listed in Table [H only (130) Elektra and (283) 
Emma have mutual orbits that are not roughly c ircular, each having an eccentricity e ~ 0.1 likely 



caused by tidal excitation ( Marchis et al. . 2008al ) 



Assuming the value of nQ of 10^^ Nm~^ for Phobos determined by Yoder ( 19821 ) is applicable 
to the components of main-belt binaries, timescales for tidal evolution are plotted in Fig. [3] and 
illustrate the interplay between material strength and system age. Systems that plot to the right 
of the 4.5 Gyr curve, the age of the Solar System, must either have smaller values than Phobos 
or have evolved very little in separation since formation. Systems that plot to the left of the 4.5 
Gyr curve are either younger than 4.5 Gyr or must have ^iQ values larger th an 10^^ N m~^ if they 
are a c tually 4.5 Gyr old. Solid rock ha s a rigidity /i of 10^'^ Nm^^ or greater ( Goldreich and Soteii . 



1966 : Dziewonski and Anderson . 198ll ) a nd, assuming Q ^ 1, the product iiQ for solid rock 
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would be of order 10^^ Nm~^ or greater; Burns and Safronov ( 19731 ) adopt a /xQ value of 3 x 10 
Nm~^ as an extreme case f or solid , non- porous rock in their analysis of damping to principal 
axis rotation. Furthermore, Harris! ( 19941 ) argues that fiQ shoul d fall within a factor of 10 of 
5 X 10^^ Nm~^^ based on analyses of Phobos and Gomet Halley ( Peale and Lissauer , 19891 ). In 
summary, a monolithic rock may have of order 10^^ N m~^ or greater while a fractured rock 
would be expected to hav e a lower txQ similar t o that of Phobos with continuing to decrease 
with increased fracturing (jHe and Ahrensl . [l993 ). 

If instead of assuming a value of /xQ for the binary systems, we set At to 4.5 Gyr, the age 
of the Solar System and the maximum amount of time these systems could have been tidally 
evolving, we can place an upper bounc^^ on jiQ with Eq. [T71 Ignoring (90) Antiope and (617) 
Patroclus, which can tidally evolve on timescales much more rapid than the age of the Solar 
System due to having nearly equal-size components, the median nQ from Table [1] f or main-belt 
binaries with 100-km-scale primaries is 1 X 10^1 Nm-2, the same value estimated bv lYodeii (jl982l ) 
for Phobos and reaso nable fo r fract ured rock. Individual values of are roughly consistent 
with the prediction of Harris! ( 1994 ) and range over more than three orders of magnitude from 
7.9 X 10* Nm~2 for (702) Alauda to 3.3 x 10^^ Nm'^ for (22) Kalliope indicating a range of 
internal structures from heavily fractured or shattered rock to solid rock if all large main-belt 
binary systems are billions of years old. If a binary system is actually younger than 4.5 Gyr, fiQ 
will scale downward by the same factor. 

Given the low fiQ values of (702) Alauda, (45) Eugenia, and (130) Elektra and the implied 
porosities of the primary components based on their densities in Table [H some of the large 



^When considering a specific value of At, Eq. [T7] provides a lower bound when ignoring the 
secondary and an upper bound when accounting for the secondary over the entire age At. Here, 
we are considering all possible values of At such that using the maximum possible age will give 
the maximum fiQ value to within a factor of two using Eq. [T7] and accounting only for tides on the 
primary component. 
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5 10 15 20 25 30 35 
Separation, o/Rp 



Fig. 3. Mass ratio q and primary-secondary separation a/Rp for 100-km-scale main-belt and Jupiter- 
Trojan binaries along with the angular momentum limit for J /J' = 0.5. The solid curve is the 
stability limit. Tidal evolution timescales are plotted assuming pp^g = 2 g cm~^, Rp = 100 km, 
and /XpQp = /^sQs = 10^^ Nm~^ and allowing tides raised on the primary only (dashed) or both 
components (dotted) to contribute over the entire evolution from an initial separation of 2Rp. With 
these parameters, (90) Antiope evolves far more rapidly than the other main-belt binaries. 
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main-belt primaries may be shattered and po rous gravitational aggregates. In the gravitational- 
aggregate model of iGoldreich and Sari void space (porosity) between the constituent 
particles, each of which has the rigidity of solid rock of order lO^'^ Nm~'^, and the increased stress 
at contact points lower the effective rigidity of the aggregate as a whole compared to monolithic 
rocl<H and can account for values of of order lO^'^ or 10^^ Nm~^ for gravitational aggregates 
100 km in radius and assuming Q ~ 10-100. Since ^jlQ values calculated for (702) Alauda, (45) 
Eugenia, and (130) Elektra fall below lO^'' Nm~'^, this could indicate that another mechanism is 
assisting tides in the semimajor-axis expansion of these systems, our ignorance of which results 
in a lower to account for the work done by the non-tidal process. If, instead, these systems 
formed near their current configurations, nQ would correspondingly increase such that the system 
evolves in semimajor axis more slowly over the same time period. However, the dependence on 
initial separation is so weak (see Section 15. 6p that the systems would have to have f ormed beyond 



lOi^p to raise //Q above 10^'' Nm ^, which is difficult by sub-catastrophic collision. iDurda et al 
find that the material excavated during a sub-catastrophic collision that initially orbits 
close to the primary typically has a semimajor axis of 4-7i2p, and further collisions among the 
orbiting material will eventually coalesce to form a single secondary orbiting within "several" 
primary radii rather than beyond 10i?p. 

The despinning timescales for the binary components of the Rs/Rp < 0.1 {q < 0.001) systems, 
assuming fiQ of 1 x 10^^ Nm~^, are of order 1 Myr for the secondaries (if they began their tidal 
evolution at an initial separation within a handful of primary radii) and much longer than the age 
of the Solar System for the primaries. Given the relatively short timescales for the secondaries, 
one would expect the secondaries in collisional binaries to have synchronized their rotation rate 
to their mean motion in the mutual orbit even if the systems are much younger than 4.5 Gyr 
old. Because the despinning timescales of the primaries far exceed the age of the Solar System, 
however, the fact the systems are still tidally evolving does not place a limit on the age of the 
binaries. The rapid synchronization of the secondary and the negligible despinning of the primary 
make our ignorance of the initial spin states of the components unimportant for tidal evolution of 
binaries with small secondaries. In addition, the small size of the secondary with respect to the 
primary and its rapid despinning timescale compared to the assumed age of the binary also make 
the correction in Eq. [T7]due to accounting for the secondary negligible [{Rg/Rp) {ts/AI) ~ 10~^ 
compared a maximum possible contribution of unity] . We also note that the smallest secondaries 
in this subs et are several kilorn eters in diameter, which gives them collisional lifetimes of more 
than 1 Gyr (|Bottke et al.l . l2005l ) such that the existence of a secondary does not limit the age of a 
large main-belt binary either. 

Taking the median fiQ of 1 x 10^^ Nm^^ found for the q < 0.1 large main-belt binaries, 
equal-mass binary (90) Antiope needs less than 10,000 years to evolve from near-contact to its 
present configuration. If the components of Antiope are heavily fractured gravitational aggregates, 
the time needed to tidally evolve becomes even shorter by an order of magnitude for every order 
of magnitude less than 10^^ Nm~^ the material strength of Antiope is. Even for monolithic rock 



Goldreich and Sari (12002) find that the minimum rigidity of a gravitational aggregate compares 



to monolithic rock as /i ~ (87rG/0^i?^/imonoiith/57eY)^^^ ~ 10^ (ii/km)Nm~^, where ey 
the yield strain. A gravitational aggregate 100 km in radius would have a minimum rigidity of 
order 10^ Nm~^. 
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with iJ,Q of order 10^'^ Nm~^, Antiope can fully evolve by tides within 1 Myr. Though Antiope is 
a perfect example of the result of tidal evolution, as it resides in a stable fully synchronous tidal 
end state, the system can evolve to its observed state so rapidly for any reasonable nQ value that 
one cannot constrain the material properties of Antiope well through tidal evolution. Timescales 
for (617) Patroclus are longer due to the wide separation of its components with tidal despinning 
acting on million-year timescales assuming the median fj,Q value. Based on its total angular 
momentum, the Patroclus system must have formed with the components separated by at least 
7Rp even if the components once rotated at the cohesionless breakup rate. However, this wider 
initial separation only affects the calculation of fxQ by 1% compared to our standard assumption 
of 2Rp. 



5.4- Near- Earth binaries 



We also consider binary systems in the near-Earth region well-characterized by radar 
and lightcurve observations, all of which have primaries with diameters o n the 1-km-sca l e or 



smaller. These system s are most likely the result of a spin -up mechanism (iMargot et al.l . 12002; 



Richardson and Walsh l » bescamns and Maxchil as evidenced by the J/ J' ~ 0.4 values 



for the near-Earth b i narie s in TableEl Since the observational confirmat ion of the asteroidal YORP 



effect (jTaylor et al.l . 120071 : iLowry et al.l . 120071: iKaasala 



nen et al.l. 120071') . spin-up via anisotropic 
200d : iBottke et al.l . bood ) has become the 



thermal re-radiation of absorbed sunlight (jRubincaml . 
preferred binary formation mechanism among km-sc ale parent bodies in the near-Earth region 
and perhaps the main belt ( Pravec and Harria . 2007). Further evidence of spin-up comes from 
the rapid rotation of the primaries that, ignoring nearly equal-mass binary (69230) Hermes, 
spin on average at ~90% of their cohesionless breakup rate with a mean period of 2.8 h and 
contain over 80% of the angular momentum of the system in their spin. Though they spin faster 
than their large main-belt counterparts, the primaries of near-Earth binaries contain a smaller 
fraction of the angular momentum of the system owing to their larger secondaries. The YORP 
spin-doubling timescale, the time required for an asteroid to change its rotation period (or rotation 
rate) by a factor of two due to radiative torques, is of order 1 Myr for a km-scale asteroid in t he 



near-Earth region (iRnbincaml . I2OO0I : IVokrouhlickv and Cape3 . 12OO2I : ICapek and Vokrouhlickyl . 



2004 ) allowing the parent body of a near-Earth binary system to spin-up rapidly and also allowing 
the post-formation primary component to maintain its rotation near the breakup rate. 

Aside from nearly equal-m ass binary (69230) Hermes that has reached its fully synchronous 
end state ( Margot et al. . 20061 ). the binaries in the near-Earth region are also in the midst of 
a lengthy tidal evolution. As shown by Fig. SI a handful of systems have evolved past the 
stability limit, but all remain far from reaching the outer synchronous orbit for J/J' = 0.4 that 
lies at 8Rp for q = 0.2 and rapidly recedes further out as q decreases. Among those beyond 
the stability limit, (35107) 1991 VH and (185851) 2000 DP107 have been well studied. (35107) 
1991 VH is believed to b e an a synchronous binary as three periods are sometimes detected in 
lightcurves ( Pravec et al.l . 20061 ). presumably the mutual orbit period, the rot ation of the p r imary , 



and the nonsynch ronous rotation of th e secondary. Tumbling of the primary (|Pravec et al.l . l200i 



or the secondary ( Margot et al. . 20081 ) could also explain the third periodicity. The orbit and 
rotation periods of the secondary of ( 185851) 2000 DPm? sugge st synchronization, but the primary 
has not been despun ( Margot et al. . 20021 : Pravec et al. . 20061 ). These systems illustrate the 
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Name 


Rp (km) 


Rs/ Rp 


Q P ( 


g cm ^) 


a/Rp 


J/J' 


(Nm-2) 


3671 


Dionysus 


0.75 


0.2 


0.0080 


2.0 


5.3 


0.35 


7.2 X 10^ 


5381 


Sekhmet 


0.5 


0.3 


0.0270 


2.0 


3.1 


0.37 


3.8 X 10^ 


35107 


1991 VH 


0.6 


0.38 


0.0549 


2.0 


6.0 


0.44 


1.4 X 10^ 


65803 


Didymos 


0.38 


0.22 


0.0106 


2.0 


3.0 


0.41 


9.8 X 10^ 


66063 


1998 ROi 


0.4 


0.48 


0.1106 


2.0 


3.6 


0.47 


3.8 X 10^ 


66391 


1999 KW4 


0.64 


0.33 


0.0359 


2.0 


4.0 


0.37 


1.6 X 10^ 


69230 


Hermes" 


0.3 


0.9 


0.7290 


2.0 


4.0 


0.49 


8.4 X 10^ 


85938 


1999 DJ4 


0.18 


0.5 


0.1250 


2.0 


4.1 


0.49 


3.7 X 10^ 


88710 


2001 SLg 


0.4 


0.28 


0.0220 


2.0 


3.8 


0.40 


5.2 X 10^ 


137170 


1999 HFi 


1.75 


0.23 


0.0122 


2.0 


3.4 


0.40 


1.1 X 10^° 


164121 


2003 YTi 


0.5 


0.18 


0.0058 


2.0 


6.4 


0.40 


6.6 X lO*' 


175706 


1996 FG3 


0.75 


0.31 


0.0298 


2.0 


3.7 


0.29 


2.7 X 10^ 


185851 


2000 DP107 


0.4 


0.41 


0.0689 


1.7 


6.6 


0.48 


2.9 X 10^ 




1994 AWi 


0.5 


0.49 


0.1176 


2.0 


4.8 


0.50 


9.7 X 10^ 




2000 UGii 


0.13 


0.58 


0.1951 


2.0 


4.3 


0.43 


2.2 X 10^ 




2002 CE26 


1.73 


0.09 


0.0007 


0.9 


2.7 


0.43 


4.4 X 10^ 




2004 DC 


0.17 


0.21 


0.0093 


1.7 


4.4 


0.41 


9.6 X 10^ 




2005 AB 


0.95 


0.24 


0.0138 


2.0 


4.0 


0.29 


1.3 X 10^ 




2005 NB7 


0.25 


0.4 


0.0640 


2.0 


3.7 


0.34 


7.5 X 10^ 




2006 GY2 


0.2 


0.2 


0.0080 


2.0 


3.0 


0.37 


2.2 X 10^ 




2007 DT103 


0.15 


0.33 


0.0359 


2.0 


3.6 


0.39 


1.6 X 10^ 



Table 2: Physical properties of near-Earth binar y asteroids with primary radii of order 1 km or 
smaher. Densities are estimated as 2 gcm~^ (see Pravec and Harrii . 2007, for details) unless the 
separation is known from radar observations, which ahows the density to be calculated via Kepler's 
Third Law. J / J' in each case roughly satisfies the condition for binary formation via spin-up of a 
single parent body. Values of /iQ correspond to tides raised on the primary only acting over the 10 
Myr dynamical lifetime of near-Earth asteroids. Accounting for tides raised on the secondary would 
increase by no more than a factor of two. Data are taken from the Ondrejov Asteroid Photom- 
etry Project binary asteroid parameters table (http:/ /www.asu.cas.cz/^asteroid/binas tdata.htm| 
2010 April 8 release) with the exceptions of 2000 DPm? from lMargot et all (j2002ll and 200413C 
whose parameters are provided by the authors. 

" (69230) Hermes is known to be in a fully synchronous state where tidal evolution has ceased. 
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necessary, but not sufficient nature of crossing the stability limit for completing tidal evolution. 
Because the inner synchronous orbit is buried within the contact limit for the g ~ 0.1 binaries 
with J I J' = 0.4 (see Fig. [T|, the secondaries have tidally evolved outward since the binaries were 
formed. 

The dynamical li fetini e of a n asteroid in the near-Earth region is of order 10 
Myr (jCladman et al.l . Il997l . l200d ). so the fact that all near-Earth binary systems lie be- 
yond the generic 10 Myr evolution curve in Fig. U] implies either the actual fiQ values are much 
weaker than the assumed 10^^ Nm~^ or that the systems must be older than 10 Myr having 
formed in the main belt prior to injection into the near-Earth region. Assuming all the binaries 
have tidally evolved for 10 Myr, the median fj.Q in Table [2] is 5 x 10^ Nm~^, more than two orders 
of magnitude smaller than the median fj,Q for the large main-belt binaries. A binary formed in 
the near-Earth region more recently than 10 Myr ago requires a smaller fiQ than presented to 
complete the same tidal evolution. If formed in the main belt, the age of the binaries could extend 
to, at most, the col l isiona l lifetime of a body 1 km in diameter, which is a few hundred million 
years (jBottke et al.l . l200,^ V Thus, the values in Table [2] could increase by roughly an order of 
magnitude if near-Earth binaries formed in the main belt. 

For the median fiQ value of 5 x 10*^ Nm~^, the secondaries in near-Earth binary systems are 
despun on 0.1-1 Myr timescales and thus are expected to be synchronously rotating for systems 
that are 10 Myr old, while the primaries are despun on timescales closer to the dynamical or 
collisional lifetimes of the bodies. Nearly equal-mass binary (69230) Hermes needs less than 1 
Myr to tidally evolve from near-contact to its fully synchronous end state assuming the median 
fj,Q value. Because the YORP timescale for km-scale asteroids in the near-Earth region is 
similar in magnitude to the tidal despinning timescale for the secondaries, YORP spin-up could 
hinder, negate, or even overcome tidal sp in-down, delaying or preventing synchronous lock of 
the secondaries ( Goldreich and Sari . I2OO9I ). The YORP timescale is much faster than the tidal 
despinning timescale for the primaries, which can help the primaries retain a rapid rotation (and 
possibly trigger a mass-lofting scenario) while tides evolve the separation of the components. 

Though far less rigid than solid rock, fi as low as 10^ N for a body 1 km in radius (10 ^ 
N m~^ for a body 100 m in radius) remain reasonable within the theory of iGoldreich and Sari 
(I2OO9I ) lor gravitational-aggregate structure (see Footnote 8). Assuming Q ~ 10-100 is applicable, 
then gravitational aggregates can account for fj,Q values of order 10^ or 10^ Nm~^ for km-scale 
near-Earth binaries. Even at such low rigidities, fi still dominates over the stress due to self-gravity 
gpR because of the low surface gravity g and small size of the typical near-Earth asteroid. 

For binary systems 2004 DC and (164121) 2003 YTi with fiQ values below lO"^ Nm'^, the 
weakness of the bodies is difficult to reconcile with Goldreich and Sari's theory of gravitational- 
aggregate structure. To justify a more rigid structure with fiQ larger than we calculate in this 
tidal evolution scenario, one of the following must be true: the system formed recently in nearly 
its current configuration leaving fiQ unconstrained, the system is older than 10 Myr, or the orbit 
expansion is aided by another mechanis m or mechan i sms, e.g., BYORP, m ass lofting, o r close 
planetary flyby. Intriguingly, 2004 DC (|Tavlor et al.l . I2OO8I ) and 2003 YTi (iNolan et al.l . I2OO4I ) 
are also both believed to have asynchronous secondaries that are not despun to the mean motion 
in their mutual orbit and have moderate eccentricities greater than 0.1. These properties place 
them in the minority among near-Earth binary systems that tend to have synchronous secondaries 
and circularized mutual orbits. Given that the secondaries are rapidly despun by tides compared 



22 




Separation, o/Rp 



Fig. 4. Mass ratio q and primary-secondary separation a/Rp for near-Earth binaries along with 
the angular momentum limit for J/ J' = 0.5. The solid curve is the stability limit. Tidal evolution 
timcscalcs are plotted assuming pp.g = 2 g cm~^, Rp = 1 km, and fJ,pQp = /isQs = 10^"*^ Nm~^ and 
allowing tides raised on the primary only (dashed) or both components (dotted) to contribute over 
the entire evolution from an initial separation of 2Rp. The spread of the binaries across a range of 
timescales indicates either a range of ages for the binaries, a range of material strengths in terms 
of fiQ, or a combination of both. 
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to orbit expansion timescales in binary systems, finding an asynchronous secondary appears to 
indicate youth, implying that the separation has not expanded much and that these binaries 
formed in essentially their current configurations. The product fiQ could be larger by an order of 
magnitude if the binaries are as old as their collisional lifetimes, but then one must explain why 
the secondaries are not despun over the longer timescale. 

The eccentric nature of the mutual orbits of 2004 DC and 2003 YTi could be key to 
understanding why these systems stand out among the rest of the nea r-Earth binary population, 
perhaps through entering a tumbling state 



saturnian satellite Hyperion (IWisdom et al 



Cuk aiid Nesvornv . 2010 simila r to that of the 



19841 ^. ICuk and Nesvorn-^ (|20inl ) suggest that after 



the secondary is initially synchronized by tides, outward BYORP evolution can both provide an 
eccentricity to the mutual orbit and trigger episodes of chaotic rotation, breaking the synchronous 
lock of the secondary. Without a synchronous secondary, BYORP must shut down in these 
systems, but BYORP may have aided tides in expanding their mutual orbits. Because of the rapid 
timescales for orbit expansion associated with BYORP, the calculated assuming only tidal 
expansion would be artificially lowered due to ignoring the additional BYORP component. In 
other words, a smaller fiQ is needed for a rapid orbital expansion, but fiQ could be larger if tides 
are not the only orbit expansion mechanism. Therefore, these asynchronous, eccentric binaries 
with very small values may be indirect evidence in favor of the BYORP scenario. While this 
scenario would help explain th e characteristics of the 20 04 DC and 2003 YTi syst e ms, on e must 
keep in mind that, contrary to Cuk and Nesvornv ( 20ld ). McMahon and Scheeres ( 2010bl ) argue 
that the BYORP effect should circularize the mutual orbit as the separation increases. Clearly, 
further research on the dynamical evolution of binary asteroid systems due to BYORP is required . 



In addition to BYORP, m ass lofting could continuously expand the orbit (jHarris et al.l . l2009i : 



Fahnestock and Scheeres . 20091 ) while a clo se planetary encounter could impulsively expand the 



orbit and provide a moderate eccentricity (jChauvineau et al.l . Il995l : IScheered - lioOTbl ). The short 
timescale for YORP spin-up could prevent some secondaries from synchronizing and may also be 
important in keeping the primary rotation rapid, hence fueling orbit expansion via mass lofting. 
If tidal despinning (and mass lofting) cannot regulate the spin of the primary, the primary could 
conceivably undergo another significant mass-shedding event produ c ing a tertiary component 
leadi ng to triple system, such as (153591) 2001 SN263 (jNolan et al.l . hooj ) and (136617) 1994 
CC (jBrozovic et al.l . I2OI0I ) in the near-Earth region, or to the ejection of one of the components 
from the system. It is likely a combination of these effects in the near-Earth region that produce 
the observed binary configurations as well as the seemingly unrealistic values of nQ for some 
systems. 



5.5. Small main-belt binaries 

Small main-belt binaries, those whose primaries have radii of a few kilometers, represent 
a population that blurs the characteristics of the large main-belt binaries and the near-Earth 
binaries. The small main-belt binaries share heliocentric orbits in the same region as their much 
larger cousins, but have similar properties to the near-Earth binaries such as angular momentum 
content. In Table [3l the average J /J' is 0.41, similar to the angular momentum required by 
the spin-fission binary formation mechanism, with 73% of the binaries having 0.3 < J/J' < 0.5. 
The minority of systems have angular momenta consistent both with sub-catastrophic collisions 
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{J / J' < 0.3) and giant impact or capture mechanisms {J/ J' > 0.5). If formed through spin-up, 
one would expect these systems to have nQ values similar to the near-Earth binaries for similar 
ages. Because these systems are discovered with lightcurve photometry, the known population is 
limited to systems with siz e ratios of roughly -Rs/-Rp > 0.2, the reliable detection limit for mutual 
occultation/eclipse events ( Pravec et al. . 20061 ). 



As seen with binaries like (90) Antiope, systems with nearly equal-mass components can 
tidally evolve to a fully synchronous end state very rapidly. If we ignore this type of system 
among the small main-belt binaries, where the spin periods have likely lengthened to equal the 
mutual orbit period, the average spin period of the primary components is 3.2 h and makes up 
~75% of the angular momentum of the system. This is similar to the 2.8-h average period for 
primaries in near-Earth binary systems and another indication that a similar binary formation 
process is at work in the near-Earth region and among km-scale main-belt asteroids. The YORP 
spin-dou bling timescale for a parent body a few km i n radius orbiting in the main belt is of order 
100 Myr (jRubincaml . I2OO0I : IVokrouhlickv and Caoekl . I2OO A ICaoek and Vokrouhlickvl . l2004h . Thus, 
in a few hundred million years, a small main-belt asteroid could spin up from a mundane period 
of 6 or 12 h to a period of 3 h or less producing a binary system with a rapidly spinning primary, 
where accounting for the time required to spin up the parent body does not drastically reduce the 
possible age of the resulting binary. The collisional lifetimes of the km-scale seconda ries, however 



limit the ages of the small main-belt binaries to between roughly 100 Myr and 1 Gyr (jBottke et al 



200,51 ). 

Assuming a Phobos-like fiQ of 10^^ Nm~^, the small main-belt binary systems in Fig. [5] 
cluster around an age of 100 Myr with an order-of-magnitude scatter to both younger and older 
ages. If small main-belt binaries have existed instead for the maximum collisional lifetime of the 
secondaries of 1 Gyr, the median fiQ in Table [3] is 7 x 10^^ Nm"^. For this value of the material 
strength, the secondaries in small main-belt binary systems are despun on 1-10 Myr timescales, 
quicker than the YORP timescale, and thus are expected to be synchronously rotating for systems 
that are 1 Gyr old, while the primaries are despun on timescales similar to the collisional lifetimes 
of the bodies or longer. Nearly equal-mass binaries can despin on 1 Myr timescales assuming the 
median nQ value. 

The median fiQ we calculate is not only larger than that found for the large main-belt 
binaries, but it is reasonable for solid or somewhat fractured rock. Considering that the majority 
of small main-belt binaries appear to have formed via spin-up, unless the primaries are essentially 
solid cores left over from a parent body that shed its loose material to form a substantial secondary 
{Rs/Rp > 0.2), the median for evolution over 1 Gyr seems like an overestimate. If the systems 
are much younger than their collisional lifetimes, for every order of magnitude younger the binary 
systems are, will also reduce by an order of magnitude. If small main-belt binaries are 100 Myr 
old, the median nQ value would be within a factor of two of Phobos or a typical large main-belt 
binary; if small main-belt binaries are 10 Myr old, they would have a median fiQ an order of 
magnitude stronger than near-Earth binaries of the same age. An order-of-magnitude difference 
in fj,Q between small main-belt binaries and near-Earth binaries could be a size effect since, in 
Goldreich and Sari's gravitational-aggregate model, the rigidity scales as the radius of the body 
(see Footnote 8) . The average small main-belt primary in Table [3] is a factor of five larger than the 
average near-Earth primary in Table [H and thus could account for the same factor in the rigidity 
of the bodies. If a size effect cannot account for the difference between the nQ values found for 
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0.40 


0.064 


2.0 


3.6 


0.36 


5.5 X 10^2 


15268 


Wendelinefroger 


2.1 


0.30 


0.027 


2.0 


5.0 


0.41 


2.8 X 10" 


16635 


1993 QO 


1.9 


0.35 


0.043 


2.0 


6.0 


0.48 


1.1 X 10" 


17260 


2000 JQss 


1.6 


0.26 


0.018 


2.0 


3.5 


0.31 


1.0 X 10^2 


26471 


2000 ASi52 


3.2 


0.36 


0.047 


2.0 


6.8 


0.42 


1.5 X 10" 


32008 


2000 HM53 


1.7 


0.50 


0.125 


2.0 


7.1 


0.50 


8.8 X 10^° 


34706 


2001 OP83 


1.6 


0.28 


0.022 


2.0 


4.4 


0.38 


2.9 X 10" 


76818 


2000 RG79 


1.4 


0.35 


0.043 


2.0 


3.5 


0.34 


2.3 X 10^2 



Table 3: Physical properties of small main-be lt binary asteroids with a, ii average primary radius of 3 
km. Densities are estimated as 2 gcm~ (see lPravec and Harrisl . I2nn7l . for details) in all cases. The 
typical J/ J' roughly satisfies the condition for binary formation via spin-up of a single parent body 
as with the near-Earth binaries. Values of correspond to tides raised on the primary acting over 
the 1 Gyr collisional lifetime of the secondaries. Accounting for tides raised on the secondary would 
increase by no more than a factor of two. Data are taken from the Ondrejov Asteroid Photom- 
etry Project binary asteroid parameters table ( ,http:/ /www.asu.cas.cz/^asteroid/binastdata.htm| 
2010 April 8 release). 



small main-belt and near-Earth binaries of similar ages, then this discrepancy could indicate that 
the //Q values found for near-Earth binaries are artificially lowered due to the presence (and our 
ignorance) of other orbit expansion mechanisms. 

Small main-belt binaries, if formed via spin-up, should form on a rather consistent basis over 
the last billion years, yet we find them to clump together at a similar age in Figure [5j For the 
median value of ^Q, this would be 1 Gyr, but the median ^jlQ value seems at odds with the weaker 
structure one expects for a binary formed via spin-up. Lowering the age of the population to 
reduce the median ^iQ of the small main-belt binaries would make it seem that the small main-belt 
binary population is missing an "older" component. For a Phobos-like fiQ of 10^^ Nm~^, binaries 
with an age of about 1 Gyr are under-represented if binaries have formed on a consistent basis, and 
the smaller the median fxQ of the binaries, the younger the observed population and also the larger 
the prop ortion of "missing" binaries. If detection bias against longer orbit periods in photometric 
surveys dPravec et al.l . hmd ) is not the sole cause of the relative dearth of wider binaries that 
would represent an older population for the same value, perhaps this is evidence that older 
(100 Myr to 1 Gyr old) binaries have been destroyed via collisions or that another mechanism 
exists in the main belt that is capable of destroying binaries. BYORP is capable of contracting 
the mutual orbits of binary systems, leading either to the coalescence of the components or a 
masking of the true age of the system by hindering tidal expansion. Unchecked orbit expansion 
by BYORP could conceivably lead to the loss of t he secondary, though con iplete dissociation of 
the components via BYORP is argued against by Cuk and Nesvorny ( 2010l ). The YORP effect 
is also capable of collapsing binaries by removing ang ular momentum from the system until a 
fully synchronous orbit no long er exists (IXavlorl . Eooi l . Continued YORP spin-up and repeated 
mass-shedding events could conceivably produce a tertiary component in less than the YORP 
spin-doubling timescale of 100 Myr. If the triple system is unstable, one body could be ejected 
such that the binaries we observe are simply the most recent incarnation of the system. Note, 
however, that while the YORP timescale is shorter than the tidal despinning timescale for the 
primary assuming the median value, if small main-belt binaries are 10 Myr old instead of 1 
Gyr old, is reduced, and tidal despinning should then overcome YORP spin- up, preventing 
further mass-shedding events and slowing down the primary's rotation from the breakup rate. 



5.6. Other considerations 



Several factors can affect the calculation of fiQ in addition to the major source of error, the 
assumed age of the binary At. These include measurement errors in the (1) sizes and densities of 
the components or (2) the current separation a/Rp, (3) the assumption of an initial separation, and 
(4) neglect of higher-order teri ns in the tidal potential at small separations. These considerations 
are discussed in more detail in lTavlor and Margo 3 (|2010l ) and are summarized here: 



The calculation of fiQ scales as Pp^'^R^ q{l + q)^^'^ oc Pp^^^ + of'^'^ for a binary with a 
known system mass such that underestimates (overestimates) in any of these parameters by 
some percentage lead to underestimates (overestimates) of /xQ by a larger percentage. The 
standard assumption of equal component densities could be troublesom e, as for (66391) 1999 
KW4, where the difference in component densities ( Ostro et al. . 20061 ) results in a 43% error 
in q that can skew pQ by a similar amount. 
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Separation, o/Rp 



Fig. 5. Mass ratio q and primary-secondary separation a/Rp for small main-belt binaries along with 
the angular momentum limit for J/ J' = 0.5. The solid curve is the stability limit. Tidal evolution 
timcscalcs arc plotted assuming pp.g = 2 g cm~^, Rp = 3 km, and fJ,pQp = /isQs = 10^^ Nm~^ and 
allowing tides raised on the primary only (dashed) or both components (dotted) to contribute over 
the entire evolution from an initial separation of 2Rp. The binaries tend to cluster near size ratios 
of 0.2 {q ~ 0.01), the lower limit for reliable mutual event detection by photometry, and an age of 
100 Myr for the assumed fiQ value. 
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2. An error in the current separation a/R^ of 10% can cause a factor of two error in the 
calculated value of //Q. Typical separation uncertainties are estimated as 10-15% by the 
Ondrejov Asteroid Photometry Project. 

3. Assuming an initial separation of 2Rp does not affect the calculated ^.Q by more than a 
factor of two unless the true initial separation is within 10% of the current (final) separation. 
We note in Tabled] that the (617) Patroclus binary system must have formed at a wider 
separation than the assumed value of 2i?p, but only if the true initial separation were greater 
than 12i?p would the estimate of /xQ increase by more than a factor of two. 

4. Allowing for higher-order terms in tidal potential due to the proximity of the components 
can increase fiQ by 20-30% for separations near 2i?p, but the contribution of extra terms 
quickly falls off to 5% for tidal evolution to 5-Rp and to 1% at lOiip. 

None of the above considerations should affect the calculation of by more than a factor of two 
or so on their own; only a severe measurement error or several of these factors working in concert 
would affect an order-of-magnitude calculation of the material properties of a binary asteroid 
system. 

Also of concern are solar tides, whose presence are felt through additional tidal distortions 
raised on each component. The ratio of the tidal amplitudes raised on the primary due to the Sun 
and the secondary scales as (Mq/Ms) (a/a©)^, where Mq is the mass of the Sun and a© is the 
heliocentric semimajor axis of the binary system. For a typical near-Earth binary, Mq/M^ ~ 10^^ 
and a/a© ~ 10"*^ such that the ratio of tidal amplitudes is 10~^ rendering solar t ides negligible 



unless the binary has a shallow perihelion of a few tenths of an astronomical unit (IScheeres et al. 



200fil : IScheeresl . l2007bl ) . Despite their small masses, the secondaries in near-Earth binary systems 



are so close to their primaries that they easily raise the dominant tide. For a typical large 
main-belt binary Mq/Ms ~ 10^^ and a/a© ~ 10~^ such that the ratio of tidal amplitudes is 10~^, 
which is still negligible in the scope of this discussion since large main-belt asteroids cannot have 
close flybys of the Sun. The wider separation of main-belt binary components weakens the tide 
raised by the secondary, but it still dominates over the solar tide. Typical small main-belt binaries 
have Mq/Ms ~ 10^^ and a/a© ~ 10~* such that the ratio of the tidal amplitudes is 10~^ making 
them the least affected by solar tides. This is due to having larger secondaries than the near-Earth 
binaries and having smaller component separations than the large main-belt binaries. Since the 
mass of the secondary is typically at least two orders of magnitude less than the primary, solar 
tides on the secondary are also negligible. 

Thus far, we have only considered binary components with equal densities. Applying a 
density ratio of 0.5 < Pp/ps < 1-5 in Eq. [7]has little effect on the locations of the fully synchronous 
orbits (see Fig. [6]) as the density ratio appears only as a pre-factor of the mass ratio q, which 
tends to be small. For typical binaries with J/ J' = 0.2 or 0.4 and q < 0.1, the effect of a density 
disparity on async/-Rp is imperceptible. For J/ J' = 0.5 and q = 1, the effect is almost 10%, similar 
to the measurement uncertainty on the separation, though one would expect similar-size bodies in 
a binary system to have similar densities. Of course, this has assumed that introducing a density 
disparity is the only alteration to Eq. [3 but changing the density ratio will necessarily change the 
mass ratio according to the definition of q (that in turn affects the calculation of pQ as in point 
7^1 above) and change the angular momentum J/ J' as well. For (66391) 1999 KW4, Pp/ps = 0.7, 
which increases the mass ratio q by 43% and increases J/ J' by 5% to support the rotation and 



29 



orbital motion of a more massive secondary. The stable fully synchronous end state shifts from 
~120i2p to ~70i?p due to the combined effects of the larger mass ratio (that shifts agync/Rp 
inward by nearly a factor of twcH) and the increased angular momentum of the system (that shifts 
a-sync/Rp outward by 10%). While this change in position is significant, the 1999 KW4 system, 
currently at 4i?p, will never attain the stable fully synchronous end state via tidal evolution due 
to the insurmountable orbit-expansion timescale at wide separations compared to the dynamical 
lifetime of a near-Earth asteroid and because the stable fully synchronous end state lies beyond 
the Hill sphere of the primary. It is clearly the effect of the density disparity on the mass ratio 
g, capable of changing g by a factor of two for 0.5 < Pp/ps < 1-5, that is more important to 
determining the locations of the stable fully synchronous tidal end state and calculating ^.Q than 
the density ratio as it appears directly in Eq. [71 

We have also only considered spherical bodies with uniform densities thus far. In Section [H 
we noted how determining the fully synchronous orbits about a nonspherical body requires a more 
complicated treatment of the mutual gravitational potential and the synchronous rotation rate 
necessary for relative equilibrium than presented here. Making a body ellipsoidal also increases 
a as it requires more angular momentum to rotate an ellipsoid at a given rate about its axis of 
maximum moment of inertia than it does a sphere; on the other hand, concentrating mass toward 
the center of the body decreases a. For a typical binary system with disparate masses, the spin of 
the primary dominates the angular momentum of the system such that a change in by some 
percentage can change the angular momentum by a similar, but lesser, percentage and affect 
the locations of the fully synchronous orbits in Fig. [1] accordingly. For systems with disparate 
masses, the spin of the secondary contributes to the total angular momentum at the 1% level or 
less such that the exact value of is unimportant. As with the density ratio, the Og/ap ratio 
only appears as a pre- factor of the mass ratio q in Eq. [71 and thus is insignificant unless g ~ 1; 
however, equal-mass binaries are dominated by the orbital rather than spin angular momentum 
when residing in a stable fully synchronous orbit. Therefore, only the value of Op is significant and 
only then for systems with components of unequal mass. Take, for example, a typical near-Earth 
binary where 80% of the angular momentum is in the spin of the primary, then allowing Qp to 
vary by 25% about the nominal value of 0.4 would in turn let the angular momentum of the 
system vary by 20% or between J/ J' = 0.32 and 0.48 assuming the nominal angular momentum 
of J / J' = 0.4 for a binary formed by spin-up. While such a range in J/ J' opens a wide swath 
in Fig. [H the stable fully synchronous orbit will typically be far from the primary, ranging from 
tens of primary radii away for g = 0.1 to hundreds or thousands of primary radii (or more) as q 
decreases, separations that are comparable to or beyond the Hill sphere of the primary. 

6. Discussion 

We have derived and plotted the locations of the fully synchronous tidal end states for all mass 
ratios and angular momenta and find that only binary asteroid systems with nearly equal-mass 
components such as (90) Antiope and (69230) Hermes reside in these end states with all other 



^In Figs. [Hand [21 for mass ratios oi q < 0.5, the location of the stable fully synchronous orbit 
rapidly recedes such that a small change in mass ratio corresponds to a drastic change in the orbital 
separation. 
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2 4 6 8 

Separation, a/Rp 



Fig. 6. Fully synchronous orbits (Eq. [7]) for binary systems with scaled angular momentum J/ J' = 
0.4 along with the contact limit (Eq. angular momentum limit (Eq. [5|), stability limit (Eq. llOp . 
and = energy limit (Eq. [T5]) for density ratios of Pp/ps = 0.5 (dashed), Pp/ps = 1 (solid), and 
Pp/Ps = 1-5 (dotted). A less dense primary shifts all curves inward, while a higher-density primary 
shifts the curves outward, though the effect is insignificant for mass ratios less than a few tenths. 
Spherical shapes with ap^s = 0.4 are assumed in all cases. 
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systems with smaller mass ratios undergoing a lengthy tidal evolution process. Though (90) 
Antiope and (69230) Hermes are excellent examples of the stable fully synchronous tidal end state, 
these systems can evolve to their current states so rapidly for reasonable values of the product of 
rigidity fx and tidal dissipation function Q that their material properties are difficult to constrain. 
Relying instead upon the binary systems with smaller mass ratios, we find among 100-km-scale 
main-belt binaries values of the material strength /xQ that, for tidal evolution over the age of the 
Solar System, are consistent with solid or fractured rock as one might expect for binaries created 
via sub-catastrophic collisions. Binaries formed in the near-Earth region via spin-up that have 
evolved over the typical 10 Myr dynamical lifetime require material strengths orders of magnitude 
smaller than their large main-belt counterparts as one might expect for gravitational aggregates of 
material compared to essentially monolithic bodies. If near-Earth binaries formed in the main belt 
prior to injection to the near-Earth region, the time available for tidal evolution can lengthen and 
fj,Q can increase in a similar manner. The discovery of many small main-belt binaries with rapidly 
spinning primaries and angular momenta similar to the near-Earth binaries lends some credence to 
the idea that near-Earth binaries, or some fraction of near-Earth binaries, could form in the main 
belt prior to delivery to the near-Earth region and have binary ages that are limited by collisional 
lifetimes rather than the dynamical lifetime in the near-Earth region. An older age for some 
near-Earth binaries, i.e., 100 Myr versus 10 Myr, would help reconcile some espec i ally s mall values 
of with the model of gravitational-aggregate structure of Goldreich and Sari ( 20091 ). If small 
main-belt binaries are up to 1 Gyr old, the maximum collisional lifetime of their secondaries, their 
material strengths are more consistent with solid rock than gravitational aggregates. Younger 
ages similar to the near-Earth binaries would reduce ^Q, but leave an older binary population 
unaccounted for either because of detection bias, collisional destruction, or another destruction 
mechanism. Smaller fiQ values for the small main-belt binaries would be more consistent with the 
gravitational-aggregate structure implied by the likely formation mechanism of YORP spin-up and 
would allow tidal spin-down to overcome continued YORP spin-up of the primary, which could 
account for the slightly slower rotation rates of small main-belt primaries compared to near-Earth 
primaries. 

Beyond the question of age, our ignorance of the inherent Q of gravitational aggregates and 
whether they can be described in such idealized terms are significant caveats in this approach 
to understanding the material properties of asteroids. Especially in the near-Earth region, the 
assumption that tides are the dominant method of mutual-orbit expansion may also be tenuous. 
The binary YORP effect, mass lofting, and close planetary encounters could conspire to also 
expand or contract the mutual orbits, either assisting or hindering the expansion of the mutual 
orbit by tides. A combination of effects altering the semimajor axis would skew the resulting 
when considering tides alone; multiple methods of orbital expansion working in concert would give 
the illusion of a small nQ value rather than the true material strength of the body. It is difficult to 
disentangle different orbit-expansion mechanisms on top of not knowing the precise age of a binary 
system. If we take tides to be the dominant orbit-expansion mechanism, then the assumption of 
the binary age At is easily the largest source of error in the calculation of fiQ for a binary system 
causing an order-of-magnitude error in fiQ for an order-of-magnitude error in At. Other sources 
of error, those summarized in Section 15.61 and also from ignoring tides raised on the secondary, 
will typically cause errors of less than a factor of two. 

In our examination of tidal evolution, a few observations have hinted at the presence of a 
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mechanism (or mec hanisms) in addition to tid es that evolves the semimajor axis of binary asteroid 
systems. FoUowing iGoldreich and Sari (|2009l '). using their model for gravitational-aggregate 
structure and assuming the actual strength of the BYORP torque on the mutual orbit is 10~^ 



times its maximum possible value i^^l it is found that BYORP should dominate the semimajor-axis 
evolution of all near-Earth binary systems with synchronous secondaries. For small main-belt 
binaries with unequal-mass components, BYORP and tides could contribute comparably to the 
semimajor-axis evolution, and even for large main-belt binaries, BYORP could be important 
under certain circumstances. Among the near-Earth binaries, the systems most susceptible to 
BYORP at some point in their evolution are 2004 DC and 2003 YTi, precisely the systems whose 
asynchronous secondaries, eccentric mutual orbits, and small and worrisome values led us to 
suspect a BYORP component to their evolution in Section 15. 4[ Large main-belt binary systems 
with mass ratios above ~0.01 are dominated by tides at their current separations. Large main-belt 
binaries with the smallest size ratios are more susceptible to BYORP, and we find that those that 
plot furthest to the right of the 4.5 Gyr curve in Fig. [3] and, hence, have the smallest calculated 

values: (45) Eugenia, (130) Elektra, and (702) Alauda, are most likely to be dominated by the 
BYORP effect at this time, which would lead to artificially low values compared to the rest of 
the binaries in Table [TJ 

In the past, at smaller separations, these low mass-ratio systems among the large main-belt 
binaries would have been tide-dominated, which would allow for synchronization of the secondary 
with a rough transition to BYORP-dominated semimajor-axis evolution occurring around QRp if 
the strength of the BYORP effect is correctly estimated. Subsequent expansion by BYORP and 
tides would then produce the currently observed systems where assuming tide-dominated evolution 
finds a smaller value for ^iQ than is necessarily true. This fact is related to point ^^3 in Section [5.61 
the calculation of is most sensitive to the evolution within 10% of the current separation such 
that if BYORP is a significant contributor to the semimajor-axis expansion in this range, then 
will be skewed to a lower value than is appropriate for the material. If the strength of the BYORP 
torque is 10~^ times its maximum possible value instead, the majority of near-Earth binaries and 
the three large main-belt systems mentioned would remain BYORP-dominated in their current 
configurations. The small main-belt binaries would be tide-dominated currently making it more 
difficult for destruction by BYORP to account for the apparent lack of older systems. If the 
strength of the BYORP torque is reduced even further to 10""^ times its maximum possible value, 
then only 2004 DC and 2003 YTi (and possibly 2000 DP107) would remain BYORP-dominated. 
Given that BYORP could have a significant impact on many of the binary systems we have 
considered, perhaps the best example for the calculation of ^iQ is (22) Kalliope, which has a small 
enough mass ratio that the system is still tidally evolving, unlike (90) Antiope or (69230) Hermes, 
yet not so small that BYORP could be an issue. Given the slew of sources of uncertainty in the 
calculation of material properties via tidal evolution, Kalliope's should lie within a factor of a 
few of 3.3 X lO^'^ Nm~^ based on a 4.5 Gyr evolution, which implies stronger material properties 
than Phobos, but is reasonable for a monolithic or fractured asteroid. 

For a tide-dominated binary system, assuming an age equal to the age of the Solar System 
provides a confident order-of-magnitude upper bound on the value of jiQ. One way to improve 



^°A similar factor for the YORP t o rque was f ound to be 4 x 10 ^ for asteroid (54509) YORP 
(formerly 2000 PH5) bv lTavlor et"all (|2007l ) and lLowrv et al.l (|2007l ). 
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on this upper bound is to limit the evolution time by using collisional lifetimes (as discussed for 
near-Earth and small main-belt binaries) or by finding binaries in collisional families. If we assume 
the binary was formed as the result of the catastrophic disruption of a p arent body that pro duced 
a family of asteroids {i.e., EEBs, escaping ejecta binaries, as described by Durda et al. ( 20041 ). then 
the binary evolves only over the age of the family. Changing the evolution time by some factor 
will affect the value of in the same manner. Of the main-belt binaries listed in Tables [J and [3l 



about a dozen (19%) are believed to be part of the collisional families deter mined bvlZappa-la et al 



(119951). smaller than the valu e of roughly 30% found among all asteroids (jMarzari et al.l . Il999l : 



Bendjoya and Zappalal. l2002l). Those families the bin aries are associated with that have age 



estirn ates, i.e., Flora (^Nesvor nv et al. . 2002al ). Themis ( Marzari et al. . 1995 . 1999 : Nesvornv et al 



20051 ) , and Eos (IVokr ouhlick y et al.l . l2006l ) , are believed to be of order 1 Gyr old so that limiting 



the ages of the binaries to the age of the families rather than the age of the Solar System does 
not change the upper bound of fj,Q by more than a factor of a few. If one could instead identify 
binary systems among mu ch younger groups, such as th e Karin or Veritas clusters, both of which 
are younger than 10 Myr (|Nesvornv et al.l . l2002bl . l2003l ). the ages of the binaries, and thus their 
fj.Q values, would be far better constrained. 
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